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ABSTRACT 

We present evidence that bona fide disks and starburst systems occupy distinct regions in the gas mass versus 
star formation (SF) rate plane, both for the integrated quantities and for the respective surface densities. This 
result is based on CO observations of galaxy populations at low and high redshifts, and on the current consensus 
for the CO luminosity to gas mass conversion factors. The data suggest the existence of two different star 
formation regimes: a long-lasting mode for disks and a more rapid mode for starbursts, the latter probably 
occurring during major mergers or in dense nuclear SF regions. Both modes are observable over a large range 
of SF rates. The detection of CO emission from distant near-IR selected galaxies reveals such bimodal behavior 
for the first time, as they allow us to probe gas in disk galaxies with much higher SF rates than are seen locally. 
The different regimes can potentially be interpreted as the effect of a top-heavy IMF in starbursts. However, we 
favor a different physical origin related to the fraction of molecular gas in dense clouds. The IR luminosity to 
gas mass ratio (i.e., the SF efficiency) appears to be inversely proportional to the dynamical (rotation) timescale. 
Only when accounting for the dynamical timescale, a universal SF law is obtained, suggesting a direct link 
between global galaxy properties and the local SF rate. 

Subject headings: galaxies: formation — cosmology: observations — infrared: galaxies — galaxies: starburst 
— galaxies: evolution 



1. INTRODUCTION 

Exploring the relation between the gas content and star 
formation rate (SFR) of galaxies is crucial to understanding 
galaxy formation and evolution. This is required to under- 
stand the nature of SF, the parameters that regulate it, and 
its possible dependence on local and global galaxy properties 
(e.g.. Silk 1997; Elmegreen 2002; Krumholz & Thompson 
2007; McKee & Ostriker 2007). In addition, this information 
is a critical ingredient of theoretical models of galaxy forma- 
tion, either based on semianalytical realizations or on numer- 
ical simulations (e.g., Guiderdoni et al. 1997; Somerville et 
al. 2001; Monaco et al. 2007; Ocvirk et al. 2008; Dekel et 
al. 2009; Gnedin et al. 2009; Croton et al. 2006). In fact, the 
physics associated with the conversion of gas into stars inside 
galaxies is overwhelmingly complicated, so that theoretical 
models generally resort to scaling laws that are calibrated us- 
in g observa tions of nearby galaxies. 

'Schmidt* (1959) first suggested the existence of a power law 
relation between surface densities of SFR and gas masses. 
Kennicutt (1998; K98 hereafter) presented a calibration of 
the Schmidt law with a slope of 1 .4 in log space, which has 
since been the most widely used in the community. K98 
fit the local populations of spiral galaxies and IR-luminous 
galaxies (LIRGs/ULIRGs), with the gas mass including both 
neutral (HI) and molecular (H2) hydrogen for spirals, and 
molecular gas only for (U)LIRGs (as their HI content is likely 
negligible). The molecular gas component is routinely esti- 
mated using its most luminous tracer, carbon monoxide (CO), 
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that is generally optically thick. This requires an empirical 
derivation of the conversion factor to derive molecular gas 
masses from CO luminosities (aco = A^gas/^-co)- 
th e Galactic value for all objects in his sample. 

iDownes & SolomonI (119981) showed that aco is smaller 
for local (U)LIRGs than for spirals, by about a factor of 
6 - see a de tailed discussion and additiona l references in 
the rev iew by lSolomon & Vanden BoutI (l2005l) . iBouche et al.l 
(l2007l) applied such different conversion factors for spirals 
and (U)LIRGs, resulting in a steeper relation with a slope 
of about 1.7. Observations of high redshift submm selected 
gal axies (SMGs) can als o be fit by the same relation, follow- 
ing iTacconi et al.l ( l2008l) who argued that a ULIRG-like aco 
is appropriate for SMGs. Integrated quantities, L'qq and IR 
luminosities (Lir ) also appea r to follow a co rrelation with a 
slope of 1 .7 (e.g.,lGreve et al.ll2005; Solo mon & Vanden Bouj 
l2005h . In the following discussion, we use a Chabrier (2003) 
IMF and a standard WMAP cosmology. 

2. RELATIONS BETWEEN IR LUMINOSITY AND GAS MASSES 

In this letter, we explore the validity of the SF law fur- 
ther by including new observations of CO emissio n in distant 
near-IR selected galaxies at z = 0.5 and z = 1.5 ( Daddi et al.l 
2008, 2010, Salmi et al. in preparation). These allow us to 
study more typical high-redshift galaxies with star formation 
rates much larger than those of local spirals but less extreme 
than those of distant SMGs. The sample of 6 CO-detected 
z = 1.5 normal (BzK-selected) galaxies is presented in DIO. 
We also use CO detections of 3 near-IR selected disk galaxies 
at z = 0.5. A detailed discussion of the z — 0.5 dataset will 
be presented elsewhere (Salmi et al. in preparation). For com- 
parison, we also show measurements for normal CO-detected 
galaxies at z = 1-2.3 from Tacconi et al. (2010), although we 
do not use these in our analysis. These new observations are 
placed in context with literature data for ULIRGs, SMGs and 
local samples of disk galaxies. 

In order to investigate the location of these populations of 
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Fig. 1. — Comparison of molecular gas masses and total IR bolometric luminosities: BzK galaxies (red filled circles; DIG), z ~ 0.5 disk galaxies (red 
filled triangles; Sal mi et al. in preparation), z = 1-2.3 normal galaxies (Tacconi et al. 2010; brown crosses), SMGs ( blue empty squ ares: Daddi et al. 2009a fj 
leaver et al. 2008; Greve et alj|2005i) . QSOs (green triangles: see Riechers et al. 200 6), local ULlRGs (black crosses ISolomon et alj ri997), local spirals (black 
filled squares, Leroy et al. 2009; black filled triangles, Wilson et al. 2009). The two nearby starbursts M82 and the nucleus of NGC 253 are also shown (data 
from Weiss et al. 2001; Houghton et al. 1997; Kaneda 2009). The solid line (Eq. 1, slope of 1.31 in the left panel) is a fit to local spirals and BzK galaxies and 
the dotted line is the same relation shifted in normalization by 1.1 dex. The dashed line in the left panel is a possible double power-law fit to spirals and BzK 
galaxies. For guidance, two vertical lines indicate SFR = 2 and 200 Mq yr^' in the right panel. 



normal high-z galaxies in the gas mass versus SFR plane, ei- 
ther for the integrated properties or for the surface densities, a 
crucial ingredient is, again, the aco conversion factor. Com- 
paring the dynamical and stellar mass estimates, DIO derive 
a high aco = 3.6 ±0.8 M© (K km s"' pc^)-'^ for the BzK 
galaxies, quite similar to that for local spirals (aco — 4.6). 
This is not unexpected, given the evidence that the z ~ 1.5 
near-IR selected galaxies appear to be high redshift analogs 
of local disks with enhanced gas content (see e.g., discus- 
sions in iDaddi et al.1 l2008l 1201 Ot iDannerbauer etaH l2009l 
Tacconi et al. 2010, and later in this letter). In the following, 
we adopt this value of aco — 3-6 for the z = 0.5-2.5 normal 
galaxies^ and the 'consensus' value for the other populations 
(aco = 4-6 for local spirals, aco = 0-8 for local (U)LIRGs 
and distant SMGs/QSOs), and explore the consequences for 
the relation between gas masses and IR luminosities/SFRs. 

Fig. [His equivalent to Fig. 13 in DIO, after replacing L^q 
with Mh2- The right panel shows the ratio of Ljr to Mm 
plotted versus Lir. The implied gas consumption timescales 
(Tgas = Mm/SFR; right panel of Fig. [D are 0.3-0.8 Gyr for 
the BzK galaxies'^, about 2-3 times that for spirals, and over 
one order of magnitude smaller for local (U)LIRGs and dis- 
tant SMGs. In a simple picture, this finding can be interpreted 
in terms of two major SF modes: a long-lasting mode ap- 
propriate for disks, that holds for both local spirals and dis- 
tant BzK galaxies, and a rapid starburst mode appropriate for 
ULIRGs, local starbursts like M82 or the nucleus of NGC253, 

' This conversion factor refers to the total gas mass, including HI, H2 and 
Helium, in their proportion within the half light radius. 
* Tacconi et al. (2010) assume a similar factor. 

' We apply a conversion SFR[My, yr"'] = 10"'" x L\^/La), treating the 
two quantities as equivalent. In the case that a significant AGN contribution 
affects Lir (e.g., for the QSOs), SFRs would be correspondingly lower. 
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Fig. 2. — The SFR density as a function of the gas (atomic and molecular) 
surface density. Red filled circles and triangles are the BzKs (DIO; filled) 
and z ~ 0.5 disks (Salmi et al. in preparation), brown crosses are z = 1- 
2.3 normal galaxies (Tacconi et al. 2010). The empty squares are SMGs: 
IBo uche et aU (2007|) (blue) and Bothwell et al. (2009) (fight green). Crosses 
and filled triangles are (U)LlRGs and spiral galaxies from the sample of K98. 
The shaded regions are THINGS spirals from Bigiel et al. (2008). The lower 
solid line is a fit to local spirals and z = 1.5 BzK galaxies (Eq. 2, slope of 
1.42), and the upper dotted fine is the same relation shifted up by 0.9 dex to 
fit local ( U)LlRGs and SM Gs. SFRs are derived from IR luminosities for the 
case of a lChabrieil <2003l) IMF. 
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Fig. 3. — The same as Fig.|2] but with the gas surface densities divided by 
the dynamical time. The best fitting relation is given in Eq. 3 and has a slope 
of 1.14. 



and distant SMGs/QSOs. 
fit: 



For the disk galaxies we formally 



log Lir/Lq = 1 .3 1 X log Mm/MQ - 2.09 



(1) 



with an error on the slope of 0.09 and a scatter of 0.22 dex. 
Combining ULIRGs and SMGs we find that they define a 
trend with a similar slope, but with about 10 times higher Ljr 
at fixed Mh2- 

A similar picture applies to the surface densities (Fig. |2]l. 
We here use the original K98 measurements for local spirals 
and (U)LIRGs, but apply our choice of aco and a Chabrier 
(2003) IMF. For consistency with the K98 relation, we mea- 
sure Egas adding HI and H2 for spirals, and H2 for IR- 
luminous galaxies, in Figs |2] and [3] The results would not 



change if we had used H2 o nly for a ll galaxies. Values for 
SMGs are taken from Bouch e et af] (|2007) 



galaxies we derive gas and SFR surface densities using the 
UV rest frame (SFR) sizes. These are consistent with the CO 
sizes (DIO) but are measured at higher S/N ratio. Again, we 
find that the populations are split in this diagram and are not 
well fit by a single sequence. Our fit to the local spirals and 
the BzK galaxies is virtually identical to the original K98 re- 
lation: 



logi:sFR/[MQyr-^kpc-^] = 1.42 x logE,,„,/ [Mq/^c-^] 



-3.83 
(2) 

the slope of 1.42 is slightly larger than that of Eq. 1, with an 
uncertainty of 0.05. The scatter along the relation is 0.33 dex. 
Local (U)LIRG and SMGs/QSOs are consistent with a re- 
lation having a similar slope and normaUzation higher by 
0.9 dex, and a scatter of 0.39 dex 

Despite their high SFR > 100 Mg yr"' and 
EsFR fZj 1A^0 yr^' kpc^^, BzK galaxies are not starbursts, 
as their SFR can be sustained over timescales comparable 
to those of local spiral disks. On the other hand, M82 



and the nucleus of NGC 253 are proto-typical starbursts, 
although they only reach a SFR of a few M0 yr^'. Following 
Figs. [T] and |2] and given the ^ 1 dex displacement of the disk 
and starburst sequences, a starburst may be quantitatively 
defined as a galaxy with Lir (or Ssfr) exceeding the value 
derived from Eqs. 1 (or 2) by more than 0.5 dex. 

The situation changes substantially when introducing 
the dynamical timescale (rrtyn) into th e picture ('Elm egreenl 
H)02; Kennicutt 1998; Kr uniiholz et alJ l2009: Silk 1991. In 
Fig. [3] we compare Sgas/^dyn to Ssfr- Measurements for 
spirals and (U)LIRGs are from K98, where r^yn is defined 
to be the rotation timescale at the galaxies' outer radius 
(although Krumholz et al. 2009 use the free fall time). For the 
near-IR/optically selected z = 0.5-2.3 galaxies we evaluate 
similar quantities at the half light radius. Extrapolating 
the measurements to the outer radius would not affect our 
results substantially. Quite strikingly, the location of normal 
high-z galaxies is hardly distinguishable from that of local 
(U)LIRGs and SMGs. All observations are well described by 
the following relation: 

\ogY.sFR/[MQyr-^kpc-^] = 

1 . 14 X \og^,ajTdyn/[MQyr-^kpc-^] - 0.62. (3) 

with a slope error of 0.03 and a scatter of 0.44 dex. The re- 
markable difference with respect to Figs.[T}{2]is due to the fact 
that the normal high-z disk galaxies have much longer dynam- 
ical timescales (given their large sizes) than local (U)LIRGs. 

We can test if this holds also for integrated quantities by 
dividing the gas masses in Fig. [T]by the average (outer ra- 
dius) dynamical timescale in each population. Spirals and 
(U)LIRGs (whose r^yn does not depend on luminosity), have 
average values of rdyn = 370 Myr and rdyn = 45 Myr, respec- 
tively (K98). To be compared to rdyn = 33 Myr for SMGs 
jBouche et al.1 120071: iTacconi et al.ll2006l) . For the QSOs, we 
use the SMG value. Assuming a flat rotation curve for BzKs, 
we get an average rdyn = 330 Myr at the outer radius, about 
3 times longer than at the half light radius, given that for an 
exponential profile 90% of the mass is enclosed within ^ 3 
half light radii. A similar value is found for our z = 0.5 disk 
galaxies and the z = 1-2.3 objects from Tacconi et al. (2010). 
Despite this simple approach, Fig.|4]shows a remarkably tight 
trend: 



For the'BzK \ogSFR/[MQyr '] = 1.42 x log(MH2/rdyn)/[M03;r 



..-ii 



0.86 
(4) 

with an error in slope of 0.05 and a scatter of 0.25 dex. Fig.|4] 
suggests that roughly 10-50% of the gas is consumed dur- 
ing each outer disk rotation for local spirals, and some 30- 
100% for BzK, z = 0.5 galaxies and local (U)LIRGs. Some 
SMGs/QSOs might even consume their gas in less than one 
rotation. 

3. DISCUSSION 

The bimodal behavior in FiglT}{2]obviously depends on our 
assumptions for aco- Even though we have assumed the 
most accurate values that are available in the literature, one 
should keep in mind that this factor is still relatively poorly 
constrained. A possible variation of aco, e.g., as a function 
of galaxy properties, may change the simple picture outlined 
here. However, based on our current understanding, this is un- 
likely to greatly affect the overall conclusions. Working only 
with observed quantities, the BzK galaxies are already offset 
from local ULIRGs in the L^q/Lir diagram by a factor of 3 
(DIO). 
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Fig. 4. — The same as Fig.[T] but dividing the gas masses by an average estimate of the outer radius rotation timescale for the different populations. The dotted 
lines in the left panel show constant fractions of gas transformed into stars in each (outer radius) rotation. The best fit in the left panel (solid line; Eq. 4) has a 
slope of 1 .42. 



In summary, the observations of near-IR selected galaxies 
at z = 1.5 do not fit a single relation of M^as versus Ljr, nei- 
ther in terms of integrated quantities nor of surface densities 
(Figs.[T]and|2l). There appears to be two main sequences, one 
for disks and one for starbursts, with the latter having 10 times 
higher Ljr at fixed Mgas (either integrated, or for surface densi- 
ties). However, all populations define a single sequence when 
dividing the gas masses or surface densities by the respective 
dynamical timescales. 

Fig.[TJ[2]cannot be interpreted in terms of a single average 
trend for all populations with a large intrinsic scatter Sig- 
nificant systematic shifts are observed between the BzK and 
z = 0.5-2.3 normal galaxies and other highly star forming 
galaxies, like local (U)LIRGs and SMGs. It is worth recalling 
here that major differences exist in the properties of local disk 
galaxies versus local (U)LIRGs, and the same is observed (at 
least on average) between distant near-IR selected BzK galax- 
ies and SMGs (see DIO, in particular Sect. 7.4, for a detailed 
discussion). Overall, this dichotomy in physical properties 
of spirals/(U)LIRGs and BzKs/SMGs makes it less surprising 
that the two classes define separate sequences in the Mgas/i-iR 
plane and correspond to two distinct SF modes. 

Previously it was argued that above 100 Mq pc^^ the 
gas conversion factor would drop by a large factor (see , 
e.g., the discussion and references in iBouche et ai] 120071: 
iTacconi et al. 2008). This does not seem to be the case for 
BzK galaxies (DIO), that exhibit SF at 1 times higher gas and 
SFR surface densities than are seen locally. One could thus 
wonder how the same SF mode can be maintained at such high 
gas surface densities. It appears that the velocity dispersion 
is high in these systems (e.g., Forster Schreiber et al. 2009), 
probably due to high turbulence. The implication would be 
that BzKs have considerably thicker disks than those of local 
spirals. A higher gas filling factor in BzK galaxies might also 
apply. 

We now discuss how the observed differences in the Mgas 
versus Ljr plane could be interpreted. We note that the di- 



chotomy could be altered in the Mgas versus SFR plane if 
stars were formed with different IMFs in disks and starbursts. 
For example, if the Lir to SFR conversion factor was lower 
for (U)LIRGs and SMGs by a factor of 10, then if we were 
to plot SFR instead of Ljr in Figs. [T]-(2] we would basically 
end up with a single SF sequence. It has been claimed that 
the IMF in b ursts could be top-heavy (e.g.. lBaugh et al.lf2005l : 
lElbaz et al.l[l995.) . While the general consensus is that the ap- 
propriate IMF for local spirals is bottom-light as in Chabrier 
(2003) or'Kroupa! (l2002h . very little is known observationally 
for the case of local (U)LIRGs and even less for high redshift 
galaxies. In principle, the IMF might explain the bimodal be- 
havior at least in part. 

On the other hand, we have shown that the ratio of Ljr to 
Mgas correlates inversely with the dynamical timescales of the 
systems. The systematic differences disappear once the gas 
masses are divided by Tdyn, given that (U)LIRGs and SMGs 
have much shorter r^yn than disk galaxies at < z < 2.3. This 
is an important result because, for the first time, we show that 
global galaxy properties (that determine rdyn) are related to 
the regulation of a local (punctual) process like SF. The dy- 
namical time is expected to scale with the gas volume density 
(p) as Tdyn °= p^"'^ (e-g-. Silk et al. 1997). This suggests that 
the different LiR/Mgas ratios are not driven by the IMF (a dif- 
ference in the IMF would break the single relations found us- 
ing Tdyn) but by the fact that in some galaxies the gas can reach 
high volume densities and the systems can have short dynam- 
ical timescales so that the fuel is consumed more rapidly and 
the resulting SFR per unit of gas mass is higher. This is con- 
sistent with the existence of a unique and approximately linear 
correlation betwe en L|r and H CN lumin osity for spirals and 
local (U)LIRGs (iGao & Solom on 2004: lJuneau et all [20091) 
and QSOs (with the exception of the highest redshift objects; 
Riechers et al. 2007; Gao et al. 2007). The HCN luminos- 
ity is a measure of dense gas, requiring H2 volume densities 
> 10"* cm^^, while CO can be thermally excited for densities 
^ 10^ cm^^. Based on the HCN/Lir relation, the bimodal 
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trends of Lir to Mgas in disks versus starbursts can be inter- 
preted in terms of a similarly bimodal behavior for the dense 
gas fractions, with roughly 10 times higher fraction of dense 
gas in the starbursts compared to disks at fixed Ljr (but also 
with about twice higher dense gas fractions in BzKs versus lo- 
cal spirals). All the observations might thus be explained by 
a genuine increase of SFR efficiency in some galaxy classes, 
probably due to the concentration of the gas at high volume 
densities. 

Major mergers or other kinds of instabilities appear to be 
the most natural explanation for the increased efficiencies in 
the starbursts (although not all mergers will necessarily pro- 
duce this effect, di Matteo et al. 2008). This is not a new 
scenario. A higher SF efficiency in merger-driven starbursts 
than in rotating disks has often been implement ed in semi- 
analytical models of galaxy f ormation (e.g., Guiderdgni et alj 
[T998; Som erville et al.ll2001l) . although this has been usually 
interpreted as the natural outcome of a single SF law with an 
exponent > 1, as long as mergers make gas lose angular mo- 
mentum and concentrate in the galaxy center. However, in 
such implementations the occurrence of very high SF rates 
in gas-rich disks is neglected; they thus have a bimodality in 
surface density, not in the SF law. Our analysis suggests that 
the original K98 calibration can account for the properties of 
disk galaxies at low and high redshifts but would underes- 
timate the SF efficiency of starbursts by a factor of 10. A 
SF law with a higher exponent of 1.7 (Bouche et al. 2007) 
would in turn overestimate SF efficiency of gas-rich disks by 
a similar amount. An implementation of such a double SF law 
would surely influence predictions from semianalytical mod- 
els of galaxy formation. 



The difference in aco for disks and starbursts helps to hide 
what we are interpreting here as large differences in the SF 
efficiency, expressed in terms of Lig^/Ma^s, by reducing the 
observed differences in Ljr/L^q. There seems to be a con- 
spiracy at work such that the particular physical conditions 
that lead to high Lm/Mgas in starbursts also determine vari- 
ations in aco that obscure observationally the differences in 
SF efficiency. We emphasize thus that the distinction of star- 
burst (or merging vs non merging) systems is important for 
interpreting CO observations, although this may be difficult 
on a case by case basis. We caution that blindly applying the 
same conversion factor to all high-z observations can lead to 
confusion. Also, care must be taken that SFR/IR-luminosities 
are accurately derived. The estimates for BzK galaxies here 
are based on the cross-comparison of 3 independent SFR in- 
dicators that agree within each other very well, and overall 
on the global assessment by Daddi et al. (2007ab) on SFR 
measuraments of near-IR selected galaxies at z ^ 2. On the 
other hand, purely radio selected or mid-IR selected popula- 
tions are likely to produce a mixed bag of merging systems 
and disk galaxies and can be affected by AGNs. 
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